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ABSTRACT 

Context. The Red MSX Source (RMS) survey is an ongoing multi-wavelength observational programme designed to return a large, 
well-selected sample of massive young stellar objects (MYSOs). We have identified ~2000 MYSOs candidates located within our Galaxy by 
comparing the colours of MSX and 2MASS point sources to those of known MYSOs. The aim of our follow-up observations is to identify 
other contaminating objects such as ultra compact (UC) HII regions, evolved stars and planetary nebulae (PNe) and distinguish between 
genuine MYSOs and nearby low-mass YSOs. 

Aims. A critical part of our follow-up programme is to conduct 13 CO molecular line observations in order to determine kinematic distances 
to all of our MYSO candidates. These distances will be used in combination with far-IR and (sub)millimetre fluxes to determine bolometric 
luminosities which will allow us to identify and remove nearby low-mass YSOs. In addition these molecular line observations will help in 
identifying evolved stars which are weak CO emitters. 

Methods. We have used the 22 m Mopra telescope, the 15 m JCMT and the 20 m Onsala telescope to conduct molecular line observations 
towards 854 MYSOs candidates located in the 3rd and 4th quadrants. These observations have been made at the 7=1-0 (Mopra and Onsala) 
and .7=2-1 (JCMT) rotational transition frequency of 13 CO molecules and have a spatial resolution of ~20"-40", a sensitivity of ^ 0. 1 K 
and a velocity resolution of ~0.2 km s -1 . 

Results. We detected 13 CO emission towards a total of 751 of the 854 RMS sources observed (~ 88%). In total 2185 emission components 
are detected above 3<r level (typically T' A > 0.3 K). Multiple emission profiles are observed towards the majority of these sources - 455 
sources (~60%) - with an average of ~4 molecular clouds along the line of sight. These multiple emission features make it difficult to 
assign a kinematic velocity to many of our sample. We have used archival CS (7=2-1) and maser velocities to resolve the component 
multiplicity towards 82 sources and have derived a criterion which is used to identify the most likely component for a further 202 multiple 
component sources. Combined with the single component detections we have obtained unambiguous kinematic velocities towards 580 sources 
(~80% of the detections). The 171 sources for which we have not been able to determine the kinematic velocity will require additional line 
data. Using the rotation curve of Brand and Blitz (1993) and their radial velocities we calculate kinematic distances for all components detected. 
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1. Introduction 

1.1. Background 

Massive stars (>8 M , >10 4 L ) are responsible for some of 
the most spectacular phenomena in astrophysics. They inject 
huge amounts of energy into the interstellar medium (ISM) in 
the form of UV radiation, molecular outflows, powerful stel- 
lar winds and eventually supernova explosions. They are also 
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responsible for depositing large amounts of enriched material 
back into the ISM. Through these processes massive stars can 
have an enormous impact on their local environment and can 
influence the evolution of their host galaxies. They are also 
thought to play an important role in regulating star formation 
through the propagation of strong shocks into their surrounding 
environment. These shocks may be responsible for triggering 
subsequent generations of star formation or disrup ting condi- 
tions necessary for star formation in nearby clouds ( lElmegreen 
1998). 
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Given the profound impact massive stars have, not only 
on their local environment, but also on a galactic scale, it 
is vital that we understand the environmental conditions and 
processes involved in their birth and the earliest stages of 
their evolution. Their rarity and short evolutionary time scales 
means that massive stars are generally located much farther 
away than low-mass star formation sites. Additionally, mas- 
sive stars are known to form exclusively in clusters making it 
difficult to identify and attribute derived quantities to individ- 
ual sources. Furthermore, they reach the main sequence whilst 
still accreting material, and heavily embedded within their na- 
tal molecular cloud, hidden away beneath many magnitudes of 
visual extinction. For these reasons, until relatively recently, 
the only catalogue of massive young stellar objects (MYSOs) 
had been limited to 3 or so serendipitously detected sources 
(He nning et al.ll 1 9841) which were mostly nearby. These diffi- 
culties have resulted in our understanding of the earliest stages 
being far less developed than in the case of low mas s star for 



mation, wh ich in broad terms is well understood (IShu et al 
1987lll993h. 



The situation has improved considerably in recent years 
with a number of studies using a variety of s election criteria 
(e.g. iMolinari et al.lfl996t I Walsh et al] 1 19971: ISridharan etaT 



20021) which have identified many new MYSOs. However, all 
of these studies are based on IRAS colours which bias them to- 
wards bright, isolated sources and tend to avoid dense clustered 
environments and the Galactic mid-plane where the majority of 
MYSOs are expected to b e found - the scale height of massive 
stars is -30' dReedl2000h . The small numbers of MYSOs iden- 
tified by these studies, and the way they have been selected, 
means that they are probably not representative of the general 
MYSO population and make statistical studies of many aspects 
of massive star formation difficult. Therefore there is an obvi- 
ous need for a large, well-selected sample with sufficient num- 
bers of sources in each luminosity bin to allow statistical stud- 
ies as a function of luminosity. The Red MSX Source (RMS) 
survey is an ongoing observational programme designed to re- 
turn just such a sample with complementary data with which to 
address these issues. 



1.2. The RMS Survey 

MYSOs are m id- and far-infrared b right with luminosities 
of 10 4 -10 5 L (IWvnn-Williamslll982h . They are likely to be 
somewhat older than the hot molecular core stage and nu- 
clear fusion has probably begun, however ongoing accretion 
in some way prevents t he ionisation of the su rroundings to 
form an HII region (see Hoare & Franco! 2007 ). They are al- 
most alway s associated wi th massive bipolar molecular out- 
flows (e.g., IWu et al.ll2004l) . another indication that major ac- 
cretion is likely to be ongoing. MYSOs are also known to pos- 
sess ionised stella r winds that are weak thermal radio sources 
(~1 mJy at 1 kp dHoareEooll) . These objects can therefore be 
roughly parameterized as near and/or mid-infrared bright, have 
luminosities consistent with young O and early B-type stars 
and are relativ ely radio continuu m quiet with some exceptions 
(e.g., S:106IR: lDrewetalJll993h . 



Lumsden et al.l (120021) compared the colours of sources 
from the MSX and 2MASS point source catalogues to those 
of known MYSOs to develop a colour selection criteria from 
which we identified approximately 2000 MYSO candidates 
spread throughout the Galaxy (\b\ <5°). Sources toward the 
Galactic centre were excluded (|/| < 10°) due to confusion and 
difficulties in calculating kinematic distances. A problem with 
a colour selected sample is that the shape of the spectral en- 
ergy distribution from an optically thick cloud is insensitive to 
the type of heating source. Therefore there are several other 
types of embedded, or dust enshrouded objects, that have simi- 
lar colours to MYSOs and contaminate our sample, such as ul- 
tra compact (UC) HII regions, evolved stars and a small number 
of planetary nebulae (PNe). 

The RMS survey is a multi-wavelength programme of 
follow-up observations designed to distinguish between gen- 
uine MYSOs and these other embedded or dusty objects 
dHoare et al. 2005 ) and to compile a database of complemen- 
tary multi-wavelength data with which to study their proper- 
tiesQ These include high resolution c m continuum observa - 
tions to identify UCHII regions and PNe dUrquhart et al.l2007l) . 
mid-infrared imaging to identify genuine point sources, ob- 
tain accurate astrometry a nd avoid excluding M YSOs located 
near UCHII region s (e. g., Mottram e t al.l 120061) . near-infrared 
spectroscopy (e.g., Clarke et al.l 2006 ) to distinguish between 
MYSOs and evolved stars. 

Another crucial part of our campaign is to obtain kinematic 
distances and luminosities towards our whole sample. These 
allow us to distinguish between nearby low and intermediate 
mass YSOs from the more massive YSOs. We have therefore 
conducted a programme of molecular line observations using 
transitions of the I3 CO molecule towards all RMS sources not 
previously observed or for which good data was not available 
(see Sect. 12. II for details). Additionally, these molecular line 
observations will allow us to eliminate a large proportion of 
evolved stars since they are not generally associated with strong 
CO emission ( 12 C O (7=1-0 and 7=2-1) typically less than 1 K; 
Loup et al.lll993l) . In this paper we present the results of our 
molecular line observations towards 854 RMS sources located 
in the southern Galactic plane (180° < I < 350°). In Sect.|2]we 
outline our source selection, observational and data reduction 
procedures. We present our results and statistical analysis in 
Sect. |3] In Sect. [4] we present a summary of the results and 
highlight our main findings. 

2. Observations and data reduction 

The 13 CO transition was chosen for our observations as it is 
generally found to be only moderately optically thick and is 
thus better able to probe the bulk parameters of gas along 
the line of sight than 12 CO. Moreover, using 13 CO avoids 
many of the problems often encountered when observing with 
12 CO, which is nearly always optically thick and can often re- 
sult in complex structures with multiple components and/or 
heavily self-absorbed spectral profiles, particularly towards the 
Galactic centre. The CS transition was considered as its high 

1 http://ast.leeds.ac.uk/RMS 
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critical density (excitation threshold ~10 4 -10 5 cm 3 ) makes it 2.2. Observation procedures 



an excellent tracer of high density molecular gas. However, CS 
is much less abundant than 13 CO and as a consequence requires 
significantly longer on-source integration time to obtain an 
equivalent signal to noise. Additionally, given the lower abun- 
dance of CS there is a the danger that many of the more distant 
sources might be missed entirely. Given its higher abundance 
and low excitation threshold 13 CO is generally much stronger 



than CS and less pr oblematic than 12 CO (c.f. IJackson et al 
2006tlWuetal1l200ll) . 



2.1. Source selection 

Of our ~2000 colour selected MYSO candidates (hereafter 
RMS catalogue) 963 are located within the 3rd and 4th quad- 
rants. The large majority of our southern RMS sources were 
previously unknown in the literature. However, approximately 
a quarter are associated with bright IRAS point sources, and 
as such, might have been previously observed as part of an- 
other programme. In order to avoid re-observing these sources 
we conducted an extensive literature search using simbaeH 
This search identified a number of sources for which archival 
data was already available, the majority of these were found 
in two large observa tional programmes: a CS survey by 
Bronfman et aL ( 19961) of colour se l ected IRAS sources, and a 
CO survey bv lWouterloot & Brand Jl989b which concentrated 
mainly on IRAS sources located in the outer Galaxy. 



Bron fman et al. (1996) conducted a programme of CS 



(7=2-1) transition o bservations towards 1427 IRAS point 
sources proposed by IWood & Churchwell dl989l) as sites of 
massive star formati on. We have 259 source s in common with 
sources observed bv lBronfman et al. (1 1996b of which they re- 
port 55 as non-dete ctions. The low number of so urces in com- 
mon is because the IWood & Churchwelll dl989b criteria pref- 
erentially selects UCHII regions rather than the earlier evolu- 
tionary phases of massive star formation. The majority of these 
were excluded from our observations, but, because CS is gener- 
ally much harder to detect than 13 CO we re-observed all of their 
non-detections (with the exception of G206.7804-01.9395 for 
which there was good CO data available). In addition to the 
their non-detections we re-observed ~25% of their detections 
to allow for consistency checks and comparison of the two trac- 
ers. 



Wouterloot & Brand! dl989h used the 12 CO (7=1-0) transi- 



tion to observe 1302 IRAS selected sources located within 10° 
of the Galactic plane and between 85° < / < 280°. We find a 
total of 76 sources in common between their observations and 
our southern MYSO candidates. Some of these were excluded 
but a number were found to possess complex profiles, or were 
found to be saturated, which led to a large number being re- 
observed. Using data in the literature we were able to reduce 
the number of observations needed by over a hundred to 856. 



http://simbad.u-strasbg.fr/simbad/sim-fid. 



We made 13 CO observations towards 856 RMS sources from a 
total of 963 located in the southern Galactic plane. The vast ma- 
jority of which were observed using the 22 m Mopra telescope 
located near Coonabarabran, New South Wales, Australia^ The 
regions of the plane not accessible from Mopra were observed 
using the 15 m James Clerk Maxwell Telescope (JCMT) lo- 
cated on Mauna Kea, Hawaifl and the 20 m Onsala telescope 
located approximately 50 km south of Gotebor^§. The 13 CO 
(7=1-0) transition was used for the Mopra and Onsala obser- 
vations, however, the instruments available on the JCMT are 
unable to operate at this frequency and therefore the higher ex- 
citation 13 CO (7=2-1) transition was used. The use of these 
two different transitions should not affect our findings since 
YSO s are expected to h ave similar intensities at these frequen- 
cies (ILittle et al.l 1 19941) . The number of sources observed by 
each telescope and the details of each telescopes' observational 
set-up are summarised in Table Q] 

All of these observations were performed in position- 
switching mode, with typical on-source integration times of 
~10 minutes for both the Mopra and Onsala observations, and 
~6 minutes for the JCMT observations (however, these integra- 
tion times were reduced for stronger sources). The on-source 
integration time was split into a number of separate scans con- 
sisting of 1 minute of on- and off-source integration. Reference 
positions were offset from source positions by 1 degree in a di- 
rection perpendicular to the Galactic plane. These were chosen 
to avoid contamination of source emission from emission in the 
reference position at a similar velocity. In some cases, partic- 
ularly towards the Galactic centre, several positions needed to 
be tried before a suitable reference position was found. 

Mopra was used to observe 818 RMS sources during three 
periods between 2002 and 2005. The receiver can be tuned to 
either single or double side-band mode. The incoming signal is 
separated into two channels using a polarisation splitter, each of 
which can be tuned separately allowing two channels to be ob- 
served simultaneously. However, the large overheads involved 
in periodically retuning the second IF between the 13 CO and 
SiO maser frequency (86.2 GHz), which was used for point- 
ing corrections, meant that for the observations made before 
2004 the second IF was set to the maser frequency and was not 
retuned back to the 13 CO frequency between pointing obser- 
vations. An upgrade of the tuning system in 2004 significantly 
improved the tuning capabilities of the receiver and made it 
feasible to use both polarisations when observing our sources 

3 Mopra is operated by the Australia Telescope National Facility, 
CSIRO and the University of New South Wales during the time of the 
observations. 

4 The James Clerk Maxwell Telescope is operated by The 
Joint Astronomy Centre on behalf of the Particle Physics and 
Astronomy Research Council of the United Kingdom, the Netherlands 
Organisation for Scientific Research, and the National Research 
Council of Canada. 

5 Onsala is operated by the Swedish National Facility for Radio 
Astronomy, Onsala Space Observatory at Chalmers University of 
Technology. 
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Fig. 1. Top panel: Galactic distribution of all southern hemisphere RMS sources located in the 3rd and 4th quadrants observed in 13 CO and presented in this paper. The colours 
correspond to whether an observation resu lted in a non-detection (red), a single component was detected (blue) or multiple components are detected (green). The integrated 
12 CO emission map of iDame et al.l (120011) (grey scale) shows the Galactic distribution of molecular material. Bottom panel: Galactic longitude-vel ocity plot showing the 
velocities of observed 13 CO components as a function of Galacti c longitude. The distribu tion of molecul ar mater i al is s hown in grey scale (again from lDame et al. I l200lh for 
comparison. The location of the spiral arms taken from model by lTaylor & Cordesl (119931) and updated by lCordesI (120041) are over-plotted in red. The regions covered by these 
observations are outlined in black (see text and Table 2 for details). 
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Table 1. Observational parameters for the 13 CO observations. 



Parameter 


Mopra 


Onsala 


JCMT 


# of sources observed 


818 


9 


27 


Rest frequency (GHz) 


110.201 


110.201 


220.3986 


Total bandwidth (MHz) 


64 


80 


267 


Vel. bandpass (km s ) 


174 


224 


360 


Number of channels 


1024 


1600 


1713 


Vel. resolution (km s~') 


0.17 


0.14 


0.21 


Beam size (") 


33 


35 


21 


Date of observations 


2002-2005 


March 2003 


2003-2004 


Integration time (mins) 


10 


10 


6 


Telescope efficiency (r]f ss ) 


0.55 


0.43 


0.69 



Table 2. Systemic spectrometer velocities used at different 
Galactic longitudes to optimise available velocity bandpass. 
These were cho sen from vi s ual in spection of the longitude- 
velocity maps of lDame et al.l (I200ll) . 

Gal. longitude Centred velocity (km s -1 ) 



180-270 


50 


271-300 


20 


301-320 


-10 


321-330 


-40 


331-350 


-75 



which significantly improved the signal-to-noise and lowered 
integration times. 

Of the remaining southern sources, 27 were observed with 
the JCMT during 2003 and 2004 as part of a long-term bad 
weather backup programme and nine were observed with the 
Onsala telescope during a run in March 2003. 

Typical system temperatures obtained for both the Mopra 
and Onsala observations are between 250-450 K depending on 
weather conditions and telescope elevation, with values typ- 
ically below 300 K. The higher frequency observations con- 
ducted with the JCMT have system temperatures between 400- 
500 K. Telescope pointing was regularly checked (~2— 3 hours) 
and found to be better than 3" for the JCMT and 10" for 
the Mopra and Onsala telescopes. To correct the measured 
antenna temperatures (7^) for atmospheric absorption, ohmic 
losses and rearward spillover and scattering, a measurement 
was made of an ambi ent load (assumed to b e at 290 K) fol- 
lowing the method of iKutner & Ulich (1981); this was done 
approximately every two hours. 

In the upper panel of Fig. Q] we present a plot of the dis- 
tribution of molecular material as a function of Galact ic lon- 
gitude and latitude as traced by the iDame et al.l (120011) 12 CO 



survey. The coloured crosses indicate the RMS positions ob- 
served as part of the 13 CO observations presented here; the 
different colours indicate whether the observation resulted in 
a non-detection, detection of a single or multiple components 
by colouring the symbols red, blue and green respectively. The 
RMS sources towards which 13 CO emission has been detected 
correlate extremely well with the distribution of the 12 CO. In 
the lower panel of Fig. Q] we present a plot of the distribution 
of 12 CO as a function of Ga lactic longitude and Vlsr (again 
taken from lDame et all200ll) . We have over-plotted the veloc- 
ities of all detected components. Additionally, we have over- 
lai d the longitude-velocity of the updated sp iral arm model 
of lTavlor& CordesI ( 1 19931) by ICordesI d2004l) in red to illus- 
trate how their velocities vary as a function of longitude in our 
Galaxy. 

The velocities over which molecular clouds can be found 
in the southern Galactic plane range from approximately 
-150 km s- 1 to +100 km s -1 . Since the total Galactic ve- 
locity range (~250 km s -1 ) is larger than the total velocity 
bandpass available to our Mopra and Onsala observations it 



0.3 
1ct rms (K) 



Fig. 2. Histograms of lcr rms obtained from 6-10 minute ob- 
servations using a bin size of 0.05 K. Over 90% of the observa- 
tions have a sensitivity of better than 0. 15 K with typical values 
of -0.1K. 



was necessary to change the central velocity as a function of 
Galactic longitude in order to optimise our velocity coverage. 
We used th e Gala ctic longitude-velocity distribution maps of 
Dame et alJ (12.0011) to centre the available velocity bandpass as 
a function of longitude to cover as much of the Galaxy as pos- 
sible. The central velocities used at various Galactic longitudes 
are tabulated in Table [2] and the region covered by each are 
outlined in black on the lower panel of Fig. Q] Comparing the 
regions covered by our observations with that of the distribu- 
tion of 12 CO shows that, although our velocity bandpass was 
limited, we have obtained an almost complete coverage of the 
southern part of the Galaxy. 



2.3. Data reduction 

The basic spectral line data reduction and processing was per- 
formed using standard procedures and the software package 
specific to the individual telescopes used: DFM (Data From 
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Mopraf] for Mopra data, the X50 package for Onsala data and 
CLASS (Continuum and Line Analysis Single-dish Software^ 
for JCMT data. As previously mentioned the observation of 
each source resulting in between six and ten separate inte- 
grations. Sky emission was subtracted from these individual 
scans which were inspected to remove poor data. The remain- 
ing scans were averaged together to produce a single spectrum 
for each source and a low-order polynomial was subtracted 
from the baseline. The resulting l<x sensitivity for each spec- 
trum was typically r^O.lK per channel (see Fig. [2j. 

The spectral line parameters were obtained by fitting 
Gaussian profiles to each component present in a given spec- 
trum using XS. Where necessary, higher-order polynomi- 
als were fitted and subtracted from the baselines before the 
Gaussians were fitted. In the majority of cases a Gaussian 
profile provided a good fit to the data, however, in a signif- 
icant number of cases, line components displayed significant 
deviation from a Gaussian profile to warrant an additional 
comment. We have adopted the classification scheme used by 
Wouterloot & Brand! dl989h for their outer Galaxy 12 CO ob- 
servations of colour selected IRAS sources. They distinguish 
seven different non-Gaussian profiles which are: flat top, asym- 
metric blue or red top, blue or red asymmetry, self-absorption, 
blue and/or red wing, blue or red shoulder, self-absorption. 
However, assigning these can be a very subjective exercise and 
care must be taken when drawing conclusions from any anal- 
ysis performed. For example, in some cases the distinction be- 
tween some of the profiles can be quite arbitrary such as be- 
tween asymmetric top and shoulder profiles. Moreover, some 
emission features can exhibit several different kinds of pro- 
file simultaneously. In Fig. [3] we present some examples of the 
emission profiles observed towards RMS sources. 

In addition to the profiles just mentioned we encountered a 
number of cases where two or more separate profiles in close 
proximity overlap with each other to form complex features, 
or where the shoulder component was particularly pronounced 
and a single Gaussian produced a poor fit to the data. In these 
cases it is not clear whether these represent the superposition 
of a number of clouds along the line of sight or supersonic mo- 
tions within a particular cloud. However, their close proxim- 
ity in velocity and resulting kinematic distances (see Sect. 13.31 ) 
would suggest that if they are not associated with the same 
cloud, they are probably part of the same molecular complex. 
In some cases it is difficult, if not impossible, to distinguish 
between a profile that appears to be self-absorbed and two 
separate, but overlapping components. We identify sources as 
blended features if either source overlaps with the other above 
the half-power level. 

Finally the data were calibrated to the corrected receiver 
temperature scale 7^, which accounts for loses due to forward 
spillover and scattering, by dividing the antenna temperatures 
by the telescope efficiency where 7£ = T^/t]{ ss , and J]f ss is the 

6 A tcl/tk graphical interface written for SPC (Spectral Line 
Reduction Package) by C. Purcell. 

7 Written by P. Bergman, Onsala Space Observatory. 

8 Part of the GILDAS (Grenoble Image and Line Data Analysis 
Software) working group software. 
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Fig. 4. Histogram of the distribution of 7^ measured towards 
RMS sources. The hard cutoff under 0.25 K on the temperature 
plot is due to the sensitivity limit. The bin size is 0.2 K. 



<-> 150 



FWHM (km a" 1 ) 

Fig. 5. Histogram showing the distribution of FWHM observed 
towards RMS sources. The distribution peaks at ~2.5 km s _1 . 
Bin size used is 0.5 km s _1 . 



forward spillover and scattering efficiency (see Table[T]for indi- 
vidual telescope efficiencies). Since the main objective of these 
observations is to obtain kinematic velocities, rather than accu- 
rate line intensities, little absolute calibration was performed. 
We therefore estimate the calibration uncertainties could be as 
much as a factor of two. 

In a small number of cases some of our spectra have 
been contaminated by emission present in the reference posi- 
tion which was not evident in the individual integrations and 
only became apparent when they were averaged together. The 
data obtained for two sources (i.e., G336. 0778+00. 0361 and 
G338.492 1+00.0683) were corrupted and will not be presented 
here. 
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Fig. 6. Histogram of integrated intensity (Jt^ dV) of all com- Fig. 8. Histogram of the number of components detected to- 
ponents detected towards RMS sources. wards RMS sources reported in this paper. 



3. Results and analysis 

3.1. Detection statistics and distributions 

We detected 13 CO emission towards 751 of the 854 RMS 
sources observed (~88%) with just over a 100 non-detections. 
Multiple emission profiles are observed towards the majority 
of these sources - 455 sources (~60%) - with an average of ~4 
components found in each spectrum. A total of 2185 emission 
components are detected above 3<x (typically T* k > 0.3 K). In 
Figs. H] and [5] we present histograms of the peak temperatures 
(TV) and FWHM of all detected components respectively. The 
hard cutoff at low T* K seen in Fig.|4]is due to the observational 
detection limit. The average temperature is ~2 K and FWHM 
is ~4 km s _1 , however, both distributions are highly skewed to 
large values and consequently more typical (i.e., median) val- 
ues are ~0.7 K and 2.5 km s _1 respectively. In Fig.[6]we present 
the distribution of integrated intensities (J Tr dV = 1 .06 dv, 
where dv is the FWHM of the line) of all of the detected com- 
ponents. 

We present the parameters obtained from Gaussian fits to 
the data and derived values in Table [3] and present plots of 
the fitted spectra in Fig. 1 1 (only available in the online ver- 
sion). The table format is as follows: the RMS names, positions 
and observations Icr rms in columns 1-4; component num- 
ber and parameters extracted from Gaussian fits (i.e., Vlsr, 
and FWHM) and profile type in columns 5-9; component in- 
tegrated intensity in column 10; in columns 11-13 we present 
the Galactic centre radius and kinematic distances determined 



from the lBrand & Blita d 19931) rotation model (see Sect. |3.3| for 
details). 

The spatial distributions of the non-detections, single and 
multiple component detections can be seen in the Galactic 
longitude-latitude plot presented in the upper panel of Fig. Q] 
In this plot the multiple component detections are clearly con- 
centrated not only towards the Galactic centre but also have a 
very narrow latitude range. The distribution of the single peak 
detections are more evenly spread in Galactic longitude but 



still show a narrow latitude range. The non-detections appear to 
be evenly distibuted in both longitude and latitude. The differ- 
ent latitude distributions can be clearly seen in the histograms 
presented in Fig. [7] The distribution of the multiple compo- 
nent sources, which constitute the majority of the detections, 
is strongly peaked around the Galactic mid-plane with nearly 
all of these sources confined within two degrees of the plane, 
and the majority within one degree. The single component de- 
tections are also strongly peaked, but have a much broader lati- 
tude distribution. The non-detections appears to be evenly dis- 
tributed in latitude; these sources are most likely evolved stars 
which have similar colours to embedded YSOs but have very 
little, if any, associated molecular material. On closer examina- 
tion of the distribution of non-detections there appears to be a 
minimum at \b\ ^ 0° possibly indicating an unfortunate line of 
sight alignment between a small number of evolved stars and 
molecular clouds towards the Galactic mid-plane. However, 
these will be identified through our near infrared spectroscopic 
observations. 



3.2. Identification of multiple components 

Although using the 13 CO transition has proved successful in 
detecting emission towards potential MYSOs one significant 
drawback is the large number of components detected towards 
the majority of our sample. A histogram of the number of com- 
ponents detected towards each of RMS source is presented 
in Fig. [8] As previously mentioned we re - observ ed a number 
of sources towards which Bronfman et al. (1996) has reported 
detecting CS emission. In sources towards which only a sin- 
gle CS and 13 CO component have been detected we find their 
velocities agree to better than 1 km s -1 . Using this velocity 
correlation we are able to identify the most likely component 
associated with an RMS source where multiple components 
are detected in 18 cases using CS (7-2-1) results found i n 
the literature (i.e., iBronfman et al. Il996t iFontani et al.lbOOSh . 
Additionally, in a number of cases the associated component 
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can easily be identified as they are far stronger than any other 
component seen in the spectrum. However, identifying the cor- 
rect component is not always so obvious, and before we can as- 
sign kinematic distances to our MYSOs candidates we need to 
find a way to resolve these multi-component sources. In bottom 
panels of Fig. [3] we present two examples of multiple compo- 
nent detections. In the left hand panel it is easy to identify the 
component associated with the RMS source, however, in the 
panel on the right two of the components are approximately 
equal in strength making it difficult to identify the one associ- 
ated with the RMS source. 

Water and methanol masers are often associated with mas- 
sive star forming regions and their velocities may help identify 
the component of interest towards some of these multiple com- 
ponent detections. Maser emission can consist of many com- 
ponents spread over a wide range of velocities, typically 10- 
15 km s -1 for methan ol and up to 70 km s _1 for water masers 
(Sri dharan et al ] l2002h . ] Despite these large velocity ranges and 
number of components, the velocity of the brightest compo- 
nent correlates extremely well with the velocity of the molec- 
ular clouds with which they are associated. In order to de- 
rive a velocity criterion with which to associate masers with 
a particular mole cular component we cro ss-matched methanol 
masers from the Pestalozzi et al. ( 20051) catalogue and water 
masers from the Arcetri ca talogue (IVald ettaro et al . 2001) with 
CS detections reported by Bronfman et alj ( 1996 ). In order to 
see how the velocities of the molecular clouds and masers var- 
ied as a function of angular separation between them we used 
three search radii (120", 60" and 15") to match masers and CS 
sources. We found 233 methanol masers and 148 water masers 
in common with the CS sources within the 120" search radius. 
In Fig. [9] we present histograms of the differences in the peak 
methanol and water maser velocities compared to the CS ve- 
locities as a function of angular separation. 

These histograms reveal a fairly robust correlation, even at 
rather large angular separations, between the observed veloci- 
ties of CS emission which gives the velocity of the molecular 
cloud and the velocities of nearby masers. These plots show 
that both distributions level off at approximately the same AV 
(~15 km s _1 ) with standard deviations of ~2.5 km s _1 and 
~3.5 km s for methanol and water masers respectively for all 
of the angular distributions plotted (source with velocity differ- 
ences greater than 200 km s _I have been excluded). It is there- 
fore possible to identify the most likely CO component in cases 
where the difference in velocity between any two CO compo- 
nents is large and the velocity difference between a particular 
CO component and the maser velocities is smaller than three 
times the standard deviation (3cr a 7.5 km s _1 and 10.5 km s _1 
for methanol and water masers respectively). 

We have used these 3cr values to identify the most likely 
star forming cloud where multiple components have been de- 
tected and either a methanol or water maser has been found 
within 120". We found 78 methanol masers and 12 water 
masers associated with 79 RMS sources and are able use this 
association criterion to resolve the multiple component ambi- 
guity towards 64 sources. Using a combination of the archival 
maser data and CS data we are able to identify the correct com- 
ponent associated with an RMS source in 82 of the 455 mul- 
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Fig. 9. Histogram illustrating the difference between the veloc- 
ities of methanol (upper panel) and water masers (lower panel) 
and the velocities of their host molecular clouds as traced by CS 
(see text for details). The bin size is 2 km s and the total num- 
ber of sources used in each plot are 233 and 148 for methanol 
and water masers respectively. The unfilled histogram includes 
all maser-CS matches within a 120" search radius, over-plotted 
in light and dark grey include all matches with in a 60" and 15" 
radius respectively. The velocity range of the water-CS sources 
have been truncated to facilitate comparison between the to dis- 
tributions shown, however, only a few sources have velocity 
differences of greater than the range plotted. 



tiple component cases. Although using these other tracers of 
massive star formation has proved useful in reducing the num- 
ber of multiple component sources by ~20%, the vast major- 
ity remain unresolved. We notice that in cases where we have 
been able to associate a component with an RMS source it is 
nearly always the component with the largest integrated inten- 
sity (J Tr dV) that is identified by maser or CS velocity. In 
fact the masers and CS lines are almost exclusively associated 
with the strongest components (in 78 of the 82 cases solved) 
and only in a few cases (4) is the emission associated with 
slightly weaker components, but never less than half that of the 
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strongest component. Therefore we can use this information to 
identify the most likely component associated with the MYSO 
by picking the strongest component as long as it is more than 
twice as strong as any other component present in the data. 

By applying this criterion to our data we are able to re- 
solve the component ambiguity towards a further 202 sources, 
bringing the total number resolved to 284 or 63% of the multi- 
component detections. Together with the single component de- 
tections we have been able to determine kinematic velocities 
towards ~80% of our detections. For sources where we are 
able to resolve the component multiplicity we indicate the most 
likely component associated with an asterisk in the compo- 
nent number column of Table [3] Additional observational data 
will be required to resolve the remaining 171 sources towards 
which multi-components are detected. A number of these may 
be resolved using data from the Methanol Multibeam (MMB) 
survey f] which is currently underway and plans to survey the 
whole Galactic plane for \b\ < 2°. For the remaining sources we 
plan to resolve their component ambiguities with a programme 
of CS, water masers and CO mapping observations. 

3.3. Kinematic distances 



We have used the rotation curve of Brand & Blita (119931) to 
calculate kinematic distances and Galactocentric radii to all 
detected components. We have assumed the distance to the 
Galactic centre is 8.5 kpc and a solar velocity of 220 km s~ ! . In 
Fig. [K)]we present a histogram of the Galactic radii of all de- 
tected components, the distribution peaks between 5 and 6 kpc, 
consistent with the location of the thick ring of material which 
surrounds the centre of the Galaxy. The vast majority of sources 
in our catalogue are located within the solar circle (~80%) and 
therefore their Vlsr results in two possible distances equally 
spaced on either side of the tangent point; these are referred to 
as the near and far distances. This is known as the kinematic 
distance ambiguity problem and makes determining accurate 
distances to objects located within the solar circle difficult. The 
only sources within the solar circle that do not suffer from this 
ambiguity are sources located at the tangent point. 

The Galactocentric radii and kinematic distances deter- 
mined from the rotation model for all sources located between 
192° < I < 350° are presented in columns 11-13 of Table [3] 
We have not presented distances for sources located near the 
Galactic anti-centre since they are very uncertain. Sources with 
components outside the solar circle have only a single distance 
as do sources located at the tangent points. Errors in the dis- 
tances have been determined by shifting the component Vlsr 
by +10 km s _1 to account for peculiar motions; these result in 
typical errors of ~1 kpc. If the near and far distances for a par- 
ticular component are within 1 kpc of each other we place it 
at the tangent point. In total 119 components were found to be 
located at or near the tangent point leaving 1959 components 
with ambiguous distances. 

There are a number of methods discussed in the liter- 
ature that can be applied to solve the dista nce ambiguity 



4 6 8 10 

Galactic Centre Radius (kpc) 



Fig. 10. Histogram showing the distribution of clouds with re- 
spected to their distance from the Galactic centre. The distribu- 
tion peaks at between 5 and 6 kpc consistent with the location 
of the 5 kpc molecular ring which surrounds the Galactic cen- 
tre. The bin size used is 1 kpc. 



19801; lArava et al.l 1200 U 120021) or radio recombination lines 
in conjunction with HI absorption toward HII regions (e.g., 



Kolp ak et al.l 120031) . and HI self-absorption and molecu- 
lar li ne emission towar d s molecular clouds and MYSOs 
(e.g.. Ijackson et al.l 120021: iBusfield et al.l 120061) . The HI self- 



such as: H2CO (formaldehyde) absorption (e.g.. lDownes et al 



http://www.jb.man.ac.uk/research/methanol 



absorption method has already been successfully used to re- 
solve di stance ambiguities to wards a subset of the RMS sam- 
ple (see IBusfield et al. 2006 for details). We will use the HI 
methods to resolve the kinematic distance ambiguity towards 
the remaining sources in a subsequent paper. 

4. Summary 

We present the results of a programme of molecular line obser- 
vations of towards 854 candidate massive young stellar objects 
(MYSOs) located in the 3rd and 4th quadrant of the Galaxy. 
These have been made using the Mopra and Onsala telescopes 
and the JCMT using the 7=1-0 and 7=2-1 rotational transi- 
tions of the 13 CO molecules. We detected emission towards 
75 1 MYSOs with multiple components being seen in 60% of 
cases. The multiple detections are concentrated within 2° of 
the Galactic mid-plane and strongly peaked within 1°. The sin- 
gle component detections are also strongly peaked but have a 
slightly broader distribution in Galactic latitude \b\ ±3°. The 
non-detections are evenly distributed over a latitude range of 
±4°, consistent with them being identified as evolved stars. 

We have used archival methanol and water maser, and CS 
(7=2-1) molecular line velocities to resolve the component 
ambiguity towards 82 MYSOs candidates where multiple com- 
ponents have been detected. Using the results of our maser and 
CS associations we have derived a criterion for selecting the 
most likely component for sources where no useful tracer is 
present. In total we have been able to identify the component 
associated with our MYSO candidates in 284 cases of the 455 



10 



J. S. Urquhart et al.: 13 CO observations of MYSO candidates 



multiple detections. Combined with the single component de- 
tections we have obtained unambiguous kinematic velocities 
towards 580 sources (~80% of the detections). The 171 sources 
for which we have not been able to determine the kinematic ve- 
locity will require additional line data. 

These observations of southern sources form part of a larger 
programme of follow-up observations of near and mid-infrared 
colour selected sample of MYSO candidates. The results of ob- 
servations made towards sources in the northern Galactic plane 
will be presented in a subsequent paper. 
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Fig. 3. Example spectra observed towards RMS sources. The spectra presented in the top four rows have been chosen to illustrate 
the different kinds of emission profiles exhibited in the data and to clarify the classification scheme used for the non-Gaussian 
profiles. Moving left to right from the top left the profiles are classified as follows: Gaussian, red wing, blue shoulder, self- 
absorbed, flat-top, red asymmetic flat-top, blended and blended feature with a main and shoulder component. The bottom two 
rows show examples of multiple components seen towards the majority of our sample. 
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Fig. 7. Histograms of the Galactic latitude distribution of observed RMS sources. The unfilled histogram in each plot consists 
of all sources observed and reported in this paper. Histograms of the latitude distributions of the multiple component detections, 
single component detections and non-detections are over-plotted (grey histogram) on the left, middle and right panels respectively. 
A bin size of 1 degree has been used for all of these histograms. 



Table 3: Sample table of source parameters. 



Field Name" 


RA 


Dec 


rms 


T l h 

Id." 


Vlsr 


T R 


FWHM 


J r R dV 


Profile type 


RGC 


Near Far 


Notes 


(J2000) 


(J2000) 


(K) 




km s 


(T) 


(km s ') 


(km s" 1 K) 




(kpc) 


(kpc) (kpc) 




G300.1615-00.0877 


12:27:08.16 


-62:49:43.68 


0.07 


1 


-39.4 


2.9 


2.4 


7.4 




7.1 


4.1 (0.72) 




G300.3412-00.2190 


12:28:36.00 


-62:58:33.24 


0.10 


1 * 


-38.5 


6.7 


2.0 


14.2 


blue shoulder 


7.1 


4.2 (0.84) 


1 










2 


34.6 


1.3 


4.2 


5.8 




10.6 


12.0 (0.99) 




G300.5320-01.7926 


12:29:04.32 


-64:33:36.72 


0.07 


1 


-0.9 


0.3 


1.0 


0.3 




8.6 


8.8 (0.91) 




G300.7027+03.0398 


12:33:44.88 


-59:45:19.44 


0.08 




















G300.7221+01.2007 


12:32:50.88 


-61:35:28.68 


0.11 


1 


-43.0 


9.6 


2.2 


22.4 


blue shoulder 


7.0 


4.1 (0.19) 




G300.9785+01.1459 


12:34:58.08 


-61:39:48.24 


0.10 


1 


-43.7 


6.4 


3.3 


22.4 


red asymmetry top 


6.9 


4.2 (0.19) 




G301.0130+01.1153 


12:35:14.64 


-61:41:45.96 


0.10 


1 


-43.1 


6.1 


3.3 


21.3 


red shoulder 


7.0 


4.2 (0.36) 




G301. 1067-01.4919 


12:34:35.52 


-64:18:14.76 


0.09 




















G301. 1175+00.9627 


12:36:02.64 


-61:51:17.28 


0.20 


1 


-41.6 


6.7 


4.7 


33.4 


red shoulder 


7.0 


4.2 (0.66) 




G301. 1726+01.0034 


12:36:31.92 


-61:49:02.64 


0.10 


1 


-40.6 


5.4 


3.8 


21.8 




7.0 


4.3 (0.78) 




G301.5309-00.8231 


12:38:48.24 


-63:39:38.16 


0.07 


1 


-5.1 


3.4 


0.7 


2.5 




8.4 


0.2 (0.93) 8.7 (0.91) 




G301. 8147+00.7808 


12:41:53.76 


-62:04:13.08 


0.07 


1 


-37.4 


9.3 


2.1 


20.7 


red wing 


7.1 


4.4 (1.11) 




G302.0213+00.2542 


12:43:31.20 


-62:36:12.96 


0.08 


1 * 


-46.1 


4.0 


2.9 


12.3 




6.9 


4.3 (0.19) 


1 










2 


-37.7 


3.1 


1.2 


3.9 


blue shoulder 


7.1 


4.5 (1.11) 




G302.4546-00.7401 


12:47:08.64 


-63:36:28.80 


0.07 


1 


30.4 


0.5 


3.1 


1.6 


blended 


10.4 


12.1 (0.97) 












2 * 


35.3 


1.3 


3.2 


4.4 


blended 


10.7 


12.6 (1.00) 


1 


G302.4867-00.0308 


12:47:31.68 


-62:53:57.48 


0.10 


1 * 


-37.1 


4.4 


3.1 


14.5 


blended (main component) 


7.1 


4.5 (1.19) 


4 










2 


-34.2 


2.6 


2.1 


5.8 


blended (red shoulder) 


7.2 


3.6(1.17) 5.6(1.34) 












3 


-29.1 


0.9 


4.0 


3.8 




7.4 


2.7 (1.50) 6.5 (1.56) 




G302.5005-00.7701 


12:47:33.12 


-63:38:19.32 


0.10 


1 


30.0 


2.6 


3.2 


8.8 


blue shoulder 


10.4 


12.1 (0.96) 




G302.6604-00.7908 


12:48:59.52 


-63:39:40.68 


0.10 


1 


23.1 


0.7 


3.4 


2.5 


blended 


9.9 


11.4 (0.92) 












2 * 


26.4 


1.3 


2.1 


2.9 


blended 


10.1 


11.8 (0.94) 


1 


G303. 5353-00.5982 


12:56:50.40 


-63:27:49.32 


0.10 


1 


18.9 


3.3 


2.3 


8.0 




9.7 


11.3 (0.90) 












2 


31.6 


0.9 


4.1 


3.9 




10.5 


12.5 (0.98) 




G303. 5990-00.6524 


12:57:25.20 


-63:30:59.40 


0.08 


1 


-2.3 


1.6 


0.4 


0.7 




8.5 


9.4 (0.87) 












2 * 


19.1 


1.8 


2.5 


4.8 




9.7 


11.3 (0.90) 


1 


G303.9973+00.2800 


13:00:41.52 


-62:34:21.72 


0.10 


1 


-45.1 


1.1 


2.3 


2.7 


blended 


6.9 


4.6 (0.87) 












2 


-42.3 


2.2 


2.4 


5.6 


blended 


6.9 


4.7 (1.06) 












3 


-28.9 


0.4 


4.4 


1.9 




7.4 


2.5 (1.47) 7.0 (1.47) 












4 


30.7 


2.4 


2.5 


6.4 




10.5 


12.5 (0.97) 




G304.57 19+00.7 163 


13:05:28.08 


-62:06:38.88 


0.10 




















G304.7700-00.5193 


13:07:49.92 


-63:19:57.72 


0.11 


1 


24.6 


1.7 


1.9 


3.4 




10.0 


12.1 (0.94) 




O305 . 1 09 / +0 1 .4684 


1 1 . f\C\ .ItC O A 

13:09:36.24 


-61:19:35. lb 


0.09 
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" Observations made with the JCMT and Onsala telescope are indicated by f and +. respectively. 

Where multiple components have been detected we identify the most likely component associated with the RMS source by an ★. 
Notes: (1) strongest component; (2) methanol maser; (3) water maser; (4) and (5) CS detection by Bronfman et al. (1996) and Fontani et al. (2005). 



Table 3: continued 
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Fig. 11. CO spectra for all RMS sources towards which emission is detected. 
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